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§1  Introduction 


Part  I  of  this  paper,  namely  Sreenath,  Oh,  Krishnaprasad,  and  Marsden  [1987],  hereafter 
denoted  [I],  studied  the  Hamiltonian  structure  and  equilibria  for  interconnected  planar  rigid  bodies, 
with  the  primary  focus  being  on  the  case  of  three  bodies  coupled  with  hinge  joints.  The 
Hamiltonian  structure  was  obtained  by  the  reduction  technique,  starting  with  the  canonical 
Hamiltonian  structure  in  material  representation  and  then  quotienting  by  the  group  of  Euclidean 
motions.  For  three  bodies,  this  Hamiltonian  structure  is  as  follows  (see  Figure  1):  the  phase  space 
is  P  =  S1  x  S1  x  R3,  parametrized  by  the  two  joint  angles  021  =:  4>  and  032  =:  y  and  three 
momenta  fi.  =  (nx,  |i2,  Ji3)  (conjugate  to  the  three  angular  variables  (0j,  02,  03)  for  the  three 
bodies)  with  the  Poisson  bracket 


rf  ffi  _  fjf _ _ dg  ^  df 

1,8  3^2  Uni  3^2  J  3$ 


fJL 

af  > 

3g 

f  dg 

3g  ^ 

"  ^3; 

ay 

U^2  " 

d[l3) 

(1.1) 


This  phase  space  is  obtained  by  first  reducing  to  center  of  mass  coordinates  and  then  getting  rid  of 
rotations  via 


P  = 


T*(S!  xS'xS1)  ,  , 

— - - - - -  =  S1  x  S1  x 


Given  a  Hamiltonian  H(<j),  y,  p.lt  p.2,  |i3),  the  evolution  equations  f  =  {f,  H}  are  equivalent  to 


3H  1 

* =  w  • 

•  -  9H  3H 

_  ay  a<t>  • 

3h 

•  =-^r* 

•  _  aH  aH 

•  a^2  ’ 

•  -  dH  9H 

¥  ”  3^3  ^2  *  . 


(1.2) 
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Figure  1 


For  this  system  the  Hamiltonian  is  shown  in  [I]  to  be 

H  =  j  ( co,  Jco  >  =  j  <  n,  J~V  >  (1.3) 

where 

(coj,  0)2,  co3)  =  (Gj,  02,  83)  = 
is  the  angular  velocity  vector  and 

II  ^>3i(<l)  +  V|/) 

J  =  ^-12^)  ^2  ^-23^V^  (1.4) 

_X3i(<})  +  \|/)  ^23^)  I3 

is  the  effective  moment  of  inertia  matrix;  the  entries  are  defined  as  follows:  let  c,  b,  e,  d  be  the 
positive  distances  shown  in  Figure  1  (i.e.,  distances  between  the  centers  of  mass  and  the  hinge 
points;  assume  here  that  the  center  of  mass  of  the  central  body  is  between  the  hinge  points)  and  let 
E”  =  mjm^mj  +  m2  +  m3),  Ij,  1^,  I3  be  the  planar  moments  of  inertia  of  each  body, 

II  =  II  +  (£12  +  £l3^c^  »  I3  =  h  +  (e23  +  e13)d2  , 
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Ar  A 

I2  =  Ij  +  (£12  "h  Ej3)b  +  (^23  £13)®  2e13be  cos  ot 

be  the  augmented  moments  of  inertia,  and 

k12(<}>)  =  (ei2  +  ei3)bc  cos  <(>  +  e13  ce  cos(<|>  +  a)  , 
^3i(a)  =  e13cd  cos  a  , 

Mv)  =  (£23  +  e13)de  cos(y  -  a)  +  e13bd  cos  y  . 


Equilibrium  solutions  are  determined  by  setting  the  time  derivatives  in  (1.2)  to  zero: 


8H  _  3H  _ 
d<|>  “  dy  ~  0  ’ 

3H  9H  3H 

3^r  =  =  3^7  =  “» 1  acons,am 


(1.5) 


To  further  simplify  the  problem  we  will  assume  that  the  center  of  mass  of  the  second  body  is 
alligned  with  the  two  hinge  points;  i.e.,  that  a  =  0  (see  Figure  1).  Then  (1.5)  is  equivalent  to  the 
system 


9H 


sin(<|>  +  y)  =  —  x  sin  <J>  , 
sin  y  =  k  sin  4>  , 


(1.6) 


where 


and 


e13(b  +  e)c  +  e12bc 
e13(b  +  e)d  +  £23de 


8i3(b  +  e)  +  E^b 
Ei3d 


(1.7) 


x  = 
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as  was  shown  in  [I]. 

It  was  also  shown  in  [I]  that  there  are  always  four  or  six  equilibria,  amongst  which  are  the 
four fundamental  equilibria 


($>  Y)  =  (0,  0),  (0,  jc),  (jc,  0),  (jc,  Ji)  .  (1.8) 

When  (k,  x)  belongs  to  the  shaded  region  in  Figure  2,  i.e.,  |x-1|<k<x  +  1,  there  are  two 
other  equilibria  determined  by 


A  1  -  (K2  +  X2)  ,  K2  —  X2  —  1 

COS  <t)  =  - X— -  and  cos  y  =  - r- -  .  (1.9) 


Correspondingly,  the  pairs  (<f>,  \jf)  lie  in  the  shaded  region  of  Figure  3. 


K 


For  example,  if  m1  =  m2  =  m3,  b  =  e  =  (id,  and  c  =  Xd,  then  k  =  X  and  x  =  3(i,  so  the 
condition  for  four  equilibria,  |x-1|<k<x  +  1,  becomes  |3|i-l|<A,<3p+l.  For  instance, 
if  p  =  1,  this  condition  is  2  <  X  <  4.  Thus  as  X.  leaves  the  range  [2, 4],  the  number  of  solutions 
drops  from  six  to  four. 
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In  [I]  it  was  shown  that  the  equilibrium  (0, 0)  representing  the  straight  stretched  out  solution 
is  stable  for  all  system  parameters.  In  §2  we  study  bifurcations  of  these  equilibria  and  we  determine 
the  eigenvalue  evolution  of  these  bifurcations  and  thereby  determine  that  the  solutions  that  are  not 
formally  stable  are  not  only  unstable,  but  are  spectrally  and  hence  exponentially  unstable,  with 
non-zero  eigenvalues  of  the  linearized  equations  on  the  real  axis.  In  §3  we  use  the  version  of  the 
Weinstein-Moser  theorem  due  to  Montaldi,  Roberts,  and  Stewart  [1987]  to  show  the  existence  of 
two  families  of  periodic  orbits  (with  symmetries)  near  the  stable  equilibrium  (<J>,  \j/)  =  (0, 0);  they 
are  shown  to  be  spectrally  stable  when  they  are  nonresonant.  Finally,  in  §4  we  show  that  the 
problem  is,  in  general,  non-integrable.  This  is  done  using  the  Melnikov  method  in  the  version 
given  by  Holmes  and  Marsden  [1983]  to  show  that  the  homoclinic  orbit  present  in  the  integrable 
case  d  =  0  leads  to  transverse  homoclinic  orbits  for  small  d  *  0.  General  conditions  for 
integrability  are  not  known  to  us. 

We  believe  that  the  periodic  solutions  found  in  §4  are  related  to  travelling  waves  in  a  long 
chain  of  n  coupled  bodies  (with  torsional  springs)  or  in  the  corresponding  continuum  limit  n  -4 
oo.  This  will  be  the  subject  of  another  investigation. 


§2  Bifurcation  of  Equilibria 


In  this  section  we  relate  the  bifurcations  of  equilibria  to  the  degeneracies  of  the  Hessian  of  the 
energy  function.  This  is  used  in  the  next  section  where  we  discuss  the  stability  indices. 

First  of  all,  one  can  see  directly  from  the  equilibrium  equations,  as  in  [I],  that  a  Hamiltonian 
pitchfork-type  bifurcation  occurs  at  each  of  the  three  unstable  fundamental  equilibrium  solutions,  as 
in  Figure  4. 


kT  T-  1 


Figure  4 
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For  example,  if  m2  =  m2  =  m3,  c  =  Xd,  and  b  =  e  =  d,  the  evolution  of  equilibria  as  X, 
ranges  from  X  <  2  to  X  >  4  is  shown  in  Figure  5. 


o~~* 


— 


G=BD 


0-^0 


•o  *  o 


^  Q~7 


o- 


Q=i 


to 


V 


Figure  5 


These  bifurcations,  which  can  all  be  seen  by  direct  calculation,  will  now  be  related  to  the 
second  variation,  or  Hessian,  of  the  Hamiltonian.  The  symplectic  leaves  of  the  phase  space  P 

are  defined  by  setting  + 1^  +  P3  =  p,  a  constant.  Equilibria  for  the  system  are  exactly  critical 
points  of  H^,  the  restriction  of  H  to  P^.  At  these  points,  the  Hessian  of  H^  is  simply  the 
restriction  of  the  second  variation  of  H  to  tangent  vectors  of  P^  at  the  equilibrium  point  in 
question.  Since  the  Hamiltonian  vector  field  restricted  to  the  leaf  has  a  zero  eigenvalue  iff  the 
Hessian  does,  it  is  a  priori  clear  that  a  bifurcation  of  equilibria  occurs  only  if  the  Hessian  along  the 
leaf  has  a  zero  eigenvalue. 

The  Hessian  is  computed  at  one  of  the  fundamental  equilibria  to  have  the  form 


82H 


r1 

o 


o 

B 


(2.1) 


as  a  5x5  matrix  with  the  variables  in  the  order  (pl5  |i2,  p3,  <{>,  \j/)  restricted  to  the  subspace 
defined  by  Sjij  +  8p2  +  §p3  =  0,  where  J  is  given  by  (1.4)  and  where 


The  Dynamics  of  Coupled  Planar  Rigid  Bodies,  II 


B  =  (positive  constant) 


cos  <)>  H —  cos(<{)  +  y) 
X 


—  cos(<j>  +  y) 

T 


1 


cos(<J>  +  y) 


—  cos  y  +  —  cos(<})  +y) 
K  X 


(2.2) 


Ignoring  the  positive  constant,  we  note  that  at  the  fundamental  equilibria. 


B(0,  0)  = 


l+1 

x 


X 


2 

X 


1  1 

- f.  — 

K  X 


B(jc,  0)  = 


1 

x 


X 


K 


B(0,  it)  = 


i-i 

X 


X 


X 


K 


B(7t,  7t)  = 


I-i 

X 


X 


X 


1  1 
K  X 


(2.3) 


Since  J  is  positive  definite,  bifurcations  of  equilibria  are  determined  by  zero  eigenvalues  of  B. 
Since 


det  B(0,  0)  =  4  +  T  +  ^  >  0  ’ 


(0, 0)  never  bifurcates.  Since 


and 


det  B(7t,  0)  = 


k-x-  1 

KX 


det  B(0, 7t)  = 


1  -  k-  x 

KX 


det  B(7t,  re)  = 


x-  k-  1 

KX 


we  can  expect  these  equilibria  to  bifurcate  at  k  =  x+1,  k=1-x,  and  k  =  x  -  1,  respectively.  As 
we  saw  above,  this  is  confirmed  by  a  direct  analysis  of  the  equilibria. 

To  analyze  the  stability  of  these  equilibria  notice  first  that  the  stretched  out  state  (0,  0)  is 
always  stable,  as  we  already  know  from  [I],  For  the  state  (7t,  k),  note  that 
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det  B(rc,  n) 


x  -  k  -  1 
KX 


trace  B(je,  n) 


2k  -  xk  -  x 

XK 


> 


so  B(jt,  n)  has 


k  >  x  -  1 
k  =  x  -  1 

K  <  X-  1 


one  negative  and  one  positive  eigenvalue, 
one  negative  and  one  zero  eigenvalue, 
two  negative  eigenvalues. 


There  are  similar  statements  for  the  equilibria  (7t,  0)  and  (0, 7t)  where  x  -  1  is  replaced  by  x  +  1 
and  1  -  x,  respectively. 


Theorem  2.1  The  equilibrium  (0,  0)  is  always  stable  and  the  other  fundamental  equilibria  are 
unstable,  and  in  fact  spectrally  unstable. 


The  proof  relies  on 


Lemma  2.2  Let  A  and  B  be  two  real  n  x  n  symmetric  invertible  matrices  with  different  numbers 
of  negative  eigenvalues.  Then  the  infinitesimally  symplectic  matrix 


- 1 

o 

w 

_  J 

■  0  I* 

i 

> 

o 

l-A  0. 

I 

1 

O 

O 

w 

1 _ 

has  at  least  one  positive  ( and  so  one  negative)  real  eigenvalue. 


Proof  of  Lemma  2.2 


det 

f 

^2n  “ 

f  0 

B" 

=  det 

ai 

-BN 

<- A 

0, 

/ 

XI , 

f 


I 


=  X2n  det 


b) 


I 


J 


(notice  that  X  *  0  since  we  assume  that  A  and  B  are  invertible) 
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U 


=  X^det 


1° 


\ 


1  + 


1 

X2 


AB 


J 


=  X2n  det  ^  ^  AB  j 


=  det  (X2  +  AB)  . 


The  lemma  follows  from  this  sublemma. 

Sublemma  2.3  Under  the  same  hypothesis  as  Lemma  2.2,  the  matrix  AB  has  at  least  one 
negative  eigenvalue. 

Proof  of  Sublemma  2.3  Since  we  assume  that  A  is  invertible.  A-1  has  the  same  inertia  index 
(number  of  negative  eigenvalues)  as  A.  Now,  consider  the  1 -parameter  family  of  symmetric 
matrices  M(t)  =  tA-1  +  (l-t)B,  0  <  t  <  1.  We  know  that  the  set  of  invertible  symmetric  matrices 
has  n  +  1  components  which  are  characterized  by  inertia  indices.  Since  M(0)  and  M(l)  have 
different  indices  and  so  are  contained  in  different  components  and  (M(t)}  is  connected,  there  must 
be  some  0  <  to  <  1  for  which  M(to)  is  not  invertible,  i.e.,  there  exists  a  non-zero  vector  v  such 
that  M(to)v  =  0,  i.e., 


(t0A-i  +  (l-to)B)v  =  0  . 


Multiplying  by  A,  we  get 


ABv 


*0 

1  “  *0 


v 


Here,  -  to/(l  -  to)  is  negative  since  0  <  to  <  1.  Therefore,  AB  has  a  negative  eigenvalue.  ■ 

Remarks  1  We  can  refine  this  lemma  to  allow  one  of  the  matrices  not  to  be  invertible.  More 
specifically,  let  us  assume  that  A  is  invertible  of  the  type  (p,  q)  and  B  is  of  type  (p\  q',  r) 
where  r  is  the  number  of  zero  eigenvalues.  Then  AB  must  have  at  least  one  negative  eigenvalue 
if  p>p'  +  r  or  q>q'  +  r.  This  then  yields  Lemma  2.2  as  before. 

2  A  generalization  of  Lemma  2.2  was  pointed  out  by  J.  Howard  [1987];  similar  criteria  for 
Krein  collisions  (Hamiltonian  Hopf  bifurcations)  would  be  of  use  in  the  case  of  three  dimensional 
coupled  rigid  bodies  (cf.  Krishnaprasad,  Grossman,  and  Marsden  [1987]). 
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(3)  Some  related  results  and  applications  are  given  in  Oh  [1987]. 

Lemma  2.3  In  the  reduced  symplectic  space  P  ,  the  Hessian  of  the  reduced  Hamiltonian  has  the 
form 

V  o' 

9 

lo  B 

where  Jf  is  a  positive  definite  2x2  matrix  and  B  is  the  matrix  (22)  in  the  canonical  coordinate 
near  the  equilibria  given  by  (<j),  \} t,  Vj,  v2),  where 

\i2  -  M-i  P-3  —  M-2 

Vi  =  “2—  »  v2  =  — —  • 

Proof  of  Lemma  2.3  This  follows  from  the  fact  that  at  the  equilibria,  d2H  are  given  by  (2.1)  in 
T*(SJ  x  S1  x  S^/S1  which  is  parametrized  by  (Pj,  p^,  p^  <j>,  \j/)  and  the  fact  that 

1—2—-  — 5 — • <t>-  vj 

are  canonical  coordinates  on  near  the  equilibria;  the  latter  is  checked  directly  from  the  bracket 

(1.1)  ■ 

Proof  of  Theorem  2.1  We  know  that  in  canonical  coordinates,  the  symplectic  structure  is  given  by 

r  o  r 
.-i  o, 

and  by  Lemma  2.3,  the  Hessian  has  the  form 

V  o' 

,0  B, 

where  the  number  of  negative  eigenvalues  of  J'  is  zero  and,  as  we  illustrated  at  B(jt,n),  B  has  at 
least  one  negative  eigenvalue.  Therefore,  by  Lemma  2.2  the  linearization  of  the  Hamiltonian  vector 
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least  one  negative  eigenvalue.  Therefore,  by  Lemma  2.2  the  linearization  of  the  Hamiltonian  vector 
field  at  the  equilibria,  namely 


f  o  r 

V 

o' 

<-l  o> 

,0 

B> 

has  at  least  one  real  eigenvalue  and  so  is  spectrally  unstable  and  so  is  non-linearly  unstable.  ■ 


Finally,  we  study  the  stability  of  the  bifurcation  branches  at  the  equilibria.  By  Theorem  2.1, 
we  know  that  the  linearization  of  the  Hamiltonian  vector  field  at  (7t,  7t)  (similarly  for  (0, 7t),  (tc, 
0))  has  at  least  one  and  so  two  real  eigenvalues.  This  came  from  observing  that 


B(jt,  it) 


u 


l  \ 

x 

T  k> 


has 


det  B(7i,  tc)  >  0  if  K  <  T  —  1  , 
det  B(tc,  7t)<0  if  k  >  t  -  1  . 


Using  Lemma  2.2,  we  get  Figure  6  for  the  positions  of  eigenvalues  of  DXH(7t,  7t)  with 
respect  to  the  parameters  (t,  k). 


k  >  x-  1 


k  <  x-  1 


Figure  6 
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By  the  spectral  property  for  the  eigenvalues  of  an  infinitesimally  symplectic  matrix,  any  small 
perturbation  of  the  DXH(7t,  n)  at  k  =  x  -  1  must  have  real  eigenvalues.  Hence  the  DXH  at  the 
bifurcated  equilibria  must  have  real  eigenvalues  at  least  near  the  bifurcation  parameter.  Thus  we 
have  the  following  theorem. 

Theorem  2.4  All  the  bifurcation  branches  from  unstable  equilibria  are  themselves  unstable. 


§3  Periodic  Solutions 


We  have  seen  that  the  straight  stretched-out  state  (<j>,  \y)  =  (0,  0)  is  stable  for  all  system 
parameters.  Correspondingly,  the  second  variation  of  the  Hamiltonian  on  the  symplectic  leaf  is 
positive  definite  at  this  point.  If  co0  denotes  the  angular  velocity  of  this  solution,  its  angular 
momentum  is 


and  its  energy  is 


3 

IL  =  Hj  +  112+^3  =  a>o  X  Jy(0’  0) 

ij=  1 


3 

Eo  =  “  <  C00(1, 1,  1),  co0J(l,  1,1  )>  =  |  ©5  £  Jjj(0,  0)  (3.1) 

y=  i 

where  J  is  the  matrix  (1.4).  The  theorem  of  Weinstein  [1973],  [1978]  and  Moser  [1976]  gives  the 
following: 

Theorem  3.1  For  any  small  e  >  0,  there  are  at  least  two  distinct  periodic  orbits  near  this 
equilibrium  on  the  energy  surface  H_1(E0  +  e)  in  the  leaf  P^. 

These  periodic  motions  in  the  reduced  symplectic  manifold  produce  quasi-periodic  motion  on 
a  torus  in  the  original  phase  space  T^S^S^S1)  by  S1  symmetry. 

Theorem  3.1  does  not  directly  tell  us  properties  of  these  periodic  orbits,  such  as  their  spatial 
structure.  We  will  obtain  such  information  in  the  case  of  a  symmetric  system  by  applying  an 
equivariant  verstion  of  the  Weinstein-Moser  Theorem  due  to  Montaldi,  Roberts,  and  Stewart 
[1987]. 

We  assume  that  the  3-body  system  is  symmetric  under  the  transformation  of  configuration 
space  given  by 


(©!,  e2,  e3)  ►-*  (e3,  e2,  in  s1  x  s1  x  s1  .  (3.2) 

This  means,  in  effect,  that  bodies  1  and  3  are  mechanically  identical.  This  assumption  gives  the 
symmetries  Jn  =  J33  and  J12  =  J23  of  the  metric  J  besides  J  being  symmetric.  The 
transformation  (3.2)  on  S1  x  S1  x  S1  induces  a  Z2-action  on  the  phase  space  T*($1  x  S1  x  S1) 
and  the  Hamiltonian  (=  the  kinetic  energy)  is  invariant  under  this  action;  this  is  our  symmetry 
assumption.  Obviously,  this  ^-action  commutes  with  the  diagonal  S1 -action  and  so  it  induces  a 
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symplectic  Z2-action  on  the  reduced  space  P^.  In  canonical  coordinates  (<{>,  \j/,  Vj,  v2)  on 
where  <}>  =  02  -  0lt  =  03  -  02,  =  (|i2  -  |ii)/2,  and  v3  =  (|i3  -  |i2)/2,  this  ^-action  can  be 

written  as 

(<>,  v,  Vj,  v2)  -4  (-  xjr,  —  —  v2,  —  v3)  .  (3.3) 

Its  fixed  manifold  is  given  by 

F=  {(021,032.Vi,v2)  |  <J)  +  \|/ =  0(mod  271^!  +v2  =  0)  .  (3.4) 

We  have  the  following  general  facts  about  fixed  manifold  under  symplectic  actions: 

Proposition  3.2  Let  a  compact  Lie  group  G  act  symplectically  on  a  symplectic  manifold  P. 
Then  each  component  of  the  fixed  point  set  Fix(G)  is  a  symplectic  submanifold  of  P. 

Proof  See,  for  instance,  Guillemin  and  Sternberg  [1984].  ■ 

Proposition  3.3  Let  H  :  P  -»  R  be  a  G-invariant  Hamiltonian  and  let  XH  be  the  associated 
Hamiltonian  vector  field.  Then  XH  is  tangent  to  each  component  of  Fix(G)  and  XH  |  Fix(G)  has 
the  Hamiltonian  H  |  Fix(G). 

Proof  See  Golubitsky  and  Stewart  [1987].  ■ 

From  these  propositions,  we  see  that  F  is  a  symplectic  submanifold  (in  fact,  one  can  easily 
check  this  without  referring  to  Proposition  3.2)  and  the  restriction  of  the  Hamiltonina  H  will  give  a 
Hamiltonian  on  F.  Notice  that  F  is  actually  diffeomorphic  to  T*SJ  which  is  two  dimensional, 
and  so  the  dynamics  on  it  is  completely  integrable.  From  the  general  fact  that  zeros  of  XH  in  F 
are  zeroes  of  XH  |  p,  w  can  see  that  XH  |  F  has  two  zeroes  and  corresponding  to  (<)>,  \\f)  =  (0,  0) 
and  (k,  n)  among  the  four  fundamental  equilibria. 

Focussing  on  the  induced  dynamics  on  F  =  T*S!,  notice  that  any  level  surfaces  of  the 
reduced  Hamiltonian  in  F  is  compact  since  those  of  the  original  Hamiltonian  H  are.  Since  we 
have  proved  that  the  induced  Hamiltonian  vector  field  XH  j  F  has  one  stable  equilibrium  and  one 
unstable  one,  the  dynamics  is  qualitatively  similar  to  the  reduced  dynamics  of  the  coupled  two  body 
case.  In  particular,  we  conclude  that  the  original  Hamiltonian  system  in  P^  has  infinitely  many 
periodic  orbits  and  at  least  two  homolinic  orbits. 

We  summarize  the  above  discussions: 
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Theorem  3.4  For  the  symmetric  coupled  planar  3-body  system,  every  symmetric  initial  condition 
gives  rise  to  a  symmetric  periodic  motion  (up  to  diagonal  action  by  S1 )  unless  the  initial  energy  is 
the  same  as  the  energy  of  the  two  equilibria  (0, 0)  and  (7t,  n).  Moreover,  the  energy  surface  of 
the  unstable  equilibrium  (rc,  7t)  contain  two  homoclinic  orbits  issuing  from  it. 

From  this  theorem  we  may  expect  that  the  slightly  unsymmetric  system  will  have  chaotic 
phenomena.  This  may  be  proved  by  an  adaptation  of  the  Melnikov  method;  see  §4  below. 

We  consider  the  dynamics  on  the  energy  surface  H-1(Eq  +  e)  for  small  e>0.  We  already 
know  from  Theorem  3.1  that  this  level  surface  contains  at  least  two  distinct  periodic  orbits.  To  get 
information  on  their  symmetry,  we  use  the  following: 

Theorem  3.5  (Equivariant  Weinstein-Moser  Theorem;  Montaldi,  Roberts,  and  Stewart  [1987].) 
Let  G  be  a  group  acting  symplectically  on  a  symplectic  maniifold  (P,  Q)  and  H  be  an  invariant 
Hamiltonian.  Let  ze  P  be  a  fixed  point  for  the  corresponding  Hamiltonian  vector  field  and 
assume: 

HI  the  Hessian  d2H(z)  is  nondegenerate;  and 

H  2  d2H(zo)  restricted  to  a  resonance  subspace  \k  is  positive  definite.  (Wk  is 
the  subspace  of  TZP  that  is  the  real  part  of  the  direct  sum  of  all  the  generalized 
eigenspaces  of  eigenvalues  of  L  =  DXH(z)  that  are  multiples  of  the  purely  imaginary 
eigenvalue  A..) 

Then  for  every  isotropy  subgroup  X  of  the  GxS1  action  on  Wk,  and  for  e>0  sufficiently 
small,  there  are  at  least 


j  dim  Fix(X,  Vk)  (3.5) 

periodic  trajectories  of  XH  with  periods  near  2tc/|  X  |  and  symmetry  group  containing  X,  on  the 
energy  surface  H  =  E0  +  e. 

To  apply  this  theorem,  we  need  information  about  the  eigenvalues  and  generalized  eigenspace 
of  the  linearization  DXH(z)  of  the  Hamiltonian  vector  field  at  the  stable  equilibrium.  The 
Hamiltonian  is  quite  complicated  and  so  it  is  tedious  to  find  the  eigenvalues  and  eigenspace  directly. 
Fortunately,  we  do  not  have  to  do  this.  Instead,  we  will  fully  exploit  the  ^-symmetry  and  solve 
using  general  facts  about  symplectic  representations  (see  Guillemin  and  Sternberg  [1984]).  We 
identify  the  tangent  space  to  P^  at  the  stable  equilibrium  with  C2  by  setting 


Zj  =  ^  +  ivj  and  Z2  =  \[r  +  iv2  . 


(3.6) 
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Then  the  induced  Z2-representation  on  this  tangent  space  is  decomposed  into  irreducible 
representations;  C2  =  C0  ©  Cj  where  C0  is  the  trivial  piece  and  C3  is  the  non-trivial  piece.  In 
fact,  C0  =  { (zj,  Z2)  |  z1  +  Z2  =  0},  which  is  the  tangent  space  to  F  at  the  equilibrium  and  Cj  = 
{(zu  Z2)  |  zl-z2  =  0). 

The  group  Z2  acts  on  Cq  trivially  and  on  Cj  by  (z,  z)  »-*  (-  z,  -  z).  Since  Hessian  of  H 
at  this  equilibrium  is  positive  definite,  all  eigenvalues  of  the  linearization  of  XH  are  imaginary  and 
come  in  pairs  ±  iXj,  ±  iX^  where  and  X2  may  be  the  same.  Since  C0  is  the  tangent  space  to 
F,  the  linearization  will  have  Cq  as  a  generalized  eigenspace  corresponding  to,  say ,  ±  iXj.  Since 
we  know  that  each  generalized  eigenspace  is  symplectic  and  pairwise  orthogonal,  Cj  will  be  the 
generalized  eigenspace  of  ±  iX^. 

We  summarize  the  above  discussions: 

Proposition  3.6  The  tangent  space  at  the  stable  equilibrium  identified  with  C2  by  (3.6)  are 
decomposed  into  irreducible  pieces  of  the  induced  representation  of  2%; 

C2  =  C0  ©  Ci 
where 

C0  =  {(zb  z2)|z1+z2  =  0}  and  Cx  =  {(zj,  z2)  |  z1  -  Z2  =  0}  .  (3.7) 


Moreover,  these  irreducible  components  correspond  to  generalized  eigenspaces  of  the  linearization 
of  XH  with  the  eigenvalues  ±  iXj,  ±  iX2,  respectively,  Xl5  X^  >  0. 

From  this  proposition  we  conclude  that  the  flows  of  the  linearization  L  on  C0  and  Cj  are 
equivalent  to 


multiplicaton  by  e271^  on  Cq  and  multiplication  by  e2"1^  on  Cj  .  (3.8) 

Thus,  we  have  ^xS1  actions  on  Cq  and  Cl5  respectively. 

Next,  we  find  the  isotropy  groups  of  these  actions  on  each  of  Cq  and  Cj: 

On  Cq  the  isotropy  group  is  Z2  x  { 1 }  and  whole  space  C0  is  the  fixed  point 
space  of  real  dimension  2. 

On  Cj  the  isotropy  group  is  {—  1 }  x  {—  1 )  and  again  the  whole  space  Cj  is  the 
fixed  point  space  of  real  dimension  2. 

Therefore,  we  have  the  following  refinement  of  Theorem  3.1: 


The  Dynamics  of  Coupled  Planar  Rigid  Bodies ,  II 


19 


Theorem  3.7  For  any  small  £  >  0 ,we  have  at  least  one  periodic  orbit  with  Trf-symmetry  and  at 
least  one  periodic  orbit  with  {-I}x{-1}  symmetry  on  the  energy  surface  H-1(E0  +  e)  in  P^. 

Remarks  1  As  we  mentioned  before,  these  periodic  orbits  give  quasi-periodic  orbits  in  the 
original  phase  space  T^S1  x  S1  x  S1).  They  have  the  pictures  in  Figure  7  in  the  stick 
representation  viewed  from  a  rotating  frame,  i.e.,  up  to  the  diagonal  S^action. 


{-1}  x  {-1 } -symmetry 


2  When  ^  #  A^  one  can  apply  an  equivariant  version  of  the  Liapunov  Center  Theorem  to 
produce  smooth  families  of  periodic  orbits  with  corresponding  symmetries  bifurcating  from  the 
stable  equilibrium. 

3  We  conjecture  that  these  perioidic  solutions  are  related  to  travelling  waves  for  many 
bodies  and  in  the  continuum  limit.  ♦ 

Finally  in  this  section,  we  examine  some  aspects  of  stability.  If  we  let  <|)t  be  the  Hamiltonian 
flow,  then  the  Floquet  operator  M(u)  of  a  periodic  orbit  u(t)  with  period  T  is  defined  by 

M(u)  =  Dc^fufO)) :  Tu(0)P  -»  Tu(0)P  .  (3.9) 

If  M(u)  has  all  eigenvalues  on  the  unit  circle,  then  u  is  called  spectrally  stable.  Note  that 
M(u)  always  has  a  generalized  eigenspace  of  dimension  at  least  2  with  eigenvalue  1 1 1  because 

u'(T)  =  D<t>r(u(0))  •  u'(0)  =  u'(0)  .  (3.10) 

Now,  let  uj  and  u2  be  the  periodic  solutions  in  Theorem  3.7  whose  periods  Tj  and  T2 
are  near  2ji/]  Aj  |  and  2n/\  A2 1,  respectively.  Then  we  have  the  following  result  about  the  spectral 
stability. 
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Theorem  3.8  If  and  X2  are  non-resonant,  then  the  two  periodic  orbits  which  were  found  in 
Theorem  3.7  are  spectrally  stable  if  £  >  0  is  small. 

Proof  Note  that  the  Roquet  operator  M(Ui)  is  close  to  exp(-  TjDXH(z))  as  £  — >  0.  Note  that  Tj 
»  2n/\  X,j  |  and  exp((-  2ti/|  |)DXh(z))  has  eigenvalue  whose  corresponding  generalized 

eigenspace  is  of  dimesion  2  and  so  has  one  simple  eigenvalue  pair  which  lies  on  the  unit  circle. 
By  the  general  rigidity  of  the  behavior  of  the  eigenvalues  of  perturbations  of  a  symplectic  matrix, 
we  conclude  that  the  eigenvalues  of  M(Ui)  stay  on  the  unit  circle  if  £  >  0  is  small.  ■ 


§4  Chaotic  Solutions 


In  this  section  we  show  that  the  dynamics  of  the  three  coupled  rigid  body  system  is  not 
integrable,  having  chaotic  solutions  of  horseshoe  type.  This  is  done  using  the  Holmes  and  Marsden 
[1983]  version  of  Melnikov's  method  to  perturb  a  homoclinic  orbit  in  a  problem  with  S1 
symmetry.  There  are  several  homoclinic  orbits  one  can  use  to  perturb;  a  pair  was  described  in 
Theorem  3.4.  Here  we  perturb  the  two  body  problem  by  adding  a  third  body  near  the  center  of 
mass  of  the  second. 

We  first  need  to  derive  an  expression  for  the  Hamiltonian  that  is  written  so  the  perturbing 
terms  are  isolated.  Refer  to  Figure  8. 


Figure  8 


Let: 

O  be  the  hinge  point  of  bodies  1  and  2; 

B  be  the  hinge  point  of  bodies  2  and  3,  and  also  the  center  of  mass  of  body  2; 
a,  b,  d  be  the  vectors  between  the  centers  of  mass  and  hinge  points  of  bodies  1,  2,  3  in 
the  reference  configuration; 

R(0)  be  the  rotation  through  angle  0; 
r  be  the  vector  from  O  to  the  system  center  of  mass; 
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r°i*  r°2*  r°3  be  vectors  from  the  system  center  of  mass  to  the  body  centers  of  mass; 

01, 02, 03  be  rotation  angles  from  a  reference  configuration  to  the  current  configuration; 
Xlt  X2,  X3  be  position  vectors  for  points  in  bodies  1,  2,  3  in  the  reference 
configuration; 

Xj,  x2,  x3  be  position  vectors  for  points  in  bodies  1, 2,  3  in  the  current  configuration. 

As  in  [I],  we  have 

Xj  =  R(Gi)Xi  +  Ti  ,  (4.1a) 

r2  =  rj  +  R(0!)a  +  R(62)b  ,  (4.1b) 

r3  =  r2  +  R(03)d  ,  (4.1c) 

mr  =  m^i  +m2r2  +  m3r3  .  (4. Id) 


We  compute  the  total  kinetic  energy  as  in  [I]  as 


j  2 

H  =  X  I  trace(coiI,cof)  +  ~ 


i=  1 


1  m2  +  m3  .  .  m3 

+  Imi - 2S-i(Ria  +  R2b)-~R3d 


+  2m2 


mi  •  x  m3  • 

—  (R,a+R2b) - -  R3d 

m  1  1  m  3 


1 

+  2m3 


mi.  .  mj  +  m3  . 

— -  (Ria  +  R2b) - - - -  R3d 

m  1  1  m  ^ 


(4.2) 


Assume  p  =  0,  without  loss  of  generality,  and  that  the  reference  configuration  is  chosen  so  a,  b, 
and  d  are  parallel.  Then  we  can  write 

H  =  Hd  =  H0  +  dH1+  0(d2)  (4.3) 

where  d  =  ||  d  ||  is  our  small  parameter, 

H0  =  —  (Ij  +  ea2)co2  +  ^  G2  +  eb2)^2  +  eab  cos  <J>  tOjO^  +  -■  , 


(4.4a) 
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and 

where 


Hj  =  y(a  cos  (<}>  +  y)  to^  +  b  cos  y 


m3(mj  +  m2)  mim3  .  ,2\  m3(mi  +  m2)  2  2 

e  =  - - -  ,  Y  =  ,  and  0(d  )  =  - - - d  ©3  . 


m 


m 


m 


We  can  rewrite  H  as  H  =  !/2  ( to,  Jdco )  where  to  =  (co1(  o^,  ®3)T, 


h  = 


i!  eabcostj)  Yad  cos(<)>  +  y)^ 

eab  cos  <j>  I2  ybd  cos  y 

yad  cos(<J>  +  y)  Y*>d  cos  y  I3 


(4.4b) 


(4.5) 


and 


0  ,  m3(m,  +  m2)  , 

Ij  =  I1+£a2  ,  I2  =  I2  +  eb2  ,  I3  =  I3  +  -M: - —  d2 


m 


Write  Jd  =  J0  +  dJj  +  0(d2)  where 


Jn  = 


(  Ii  eab  cos  <{>  O') 

eab  cos  <J>  I2  0 

V  0  0  I3, 


and 


Ji  = 


0  0  ya  cos(<f>  +  y)') 

0  0  Y^  cos  y 

Vya  cos(<))  +  y)  yb  cos  y  0 


We  need  to  write  the  kinetic  energy  with  respect  to  the  momentum  jx  rather  than  the  angular 
velocity.  This  is  done  using  n  =  Jdco: 

Hd(|i)  =  =  H0(^)  +  dH1(^)  +  O(d2) 
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where 


3u  "ji 

Hi<W  =  St  =  -<M»'  Si  a.o10' V>  =  -<!*•  V'JilrtD 


(4.6) 


Here, 


Jo”1 


f  f 
A 


I-j 


eab  cos  <J> 
0 


-  eab  cos 

h 

0 


I 


0  ^ 
0 

1 

3  J 


(where  A  =  —  c  a  b  cos  0),  and  so 


Jo^W  = 


0 


0 


0 


0 


*3 


-l 


ft''  „  _  2. 


r-1 


(?2a  C°s(<t)  +  V) 
-  eb  a  cos  <J>  cos  \j /) 


r-l 


yl3  2 

— [(-  ea  b  cos(<{)  +  y) 

•  cos  <>)  +  I3b  cos  y)] 


(I2a  cos(<{)  +  y)  — [(-  ea  b  cos(<))  +  y) 
eb2a  cos  <{>  cos  y)  •  cos  <}> )  + 1 1  b  cos  y)3 


0 


Therefore, 


H°  =  2 J  ^ 1  +  W“  ^  C0S  ^  +  \  13  W 


(4.7a) 


and 
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yl3  *a  2 

Hj  = - ^ cos(4>  +  y)  -  eb  cos  <])  cos  y) 


Ylj-'b 


|X2ji3(-  ea  cos(<J)  +  y)  cos  <{>  +  Ij  cos  y)  . 


(4.7b) 


Now,  notice  that  when  d  =  0,  i.e.,  B  coincides  with  the  center  of  mass  of  body  3,  the  system  is 
completely  integrable  and  we  know  that  the  reduced  system  has  two  homoclinic  orbits,  given  in 
Figure  1  of  [I].  This  system  is  in  the  framework  of  the  Melnikov  method  with  an  S1  symmetry  as 
generalized  by  Holmes  and  Marsden  [1983]. 

From  (4.7a)  we  see  that  the  unperturbed  Hamiltonian  Hq  has  an  additional  S1  symmetry 
given  by  y-rotations;  or,  equivalently,  by  03-rotations  in  the  original  system,  which  induces  the 
obvious  Poisson  action  whose  moment  mapping  is  exactly  p3  :  T*(S!  x  S1  x  S^/S1  —>  R.  Note 
that  Pi  +  p2  +  p3  is  the  moment  mapping  corresponding  to  the  simultaneous  uniform  rotation;  i.e., 
the  diagonal  action  of  S1.  The  Hd  and  Ho-flows  restrict  to  the  symplectic  leaf  PM  =  {pl5  p2,  P3, 
<J),  y  |  Pj  +p2  +  p3  =  M),  and  (y,  p3)  are  conjugate  variables  in  this  symplectic  leaf.  Note  that  the 
equations  of  motion  for  Ho  are  given  by 


Pi  = 


9H0 

9H0 

3H0 

90 

>  Ms 

9<J> 

+  ’ 

M3  = 

<j>  = 

9H0 

3H„ 

.  3H0 

9Ho 

^2 

9pi 

* 

V=  3S7 

^2 

dHo 

9y 


(4.8) 


By  regrouping  equations  (4.8),  we  see  that  (p1?  p2,  <J>)  can  be  separated;  after  solving  this  system 
we  can  substitute  back  to  get  the  equations  for  (p3,  y).  Since  y  is  the  cyclic  variable  for  Hq, 


p3  =  0  ,  y  =  I3_1p3  -  ~  (I1P2  “  eab  cos  <(>  Pj)  (=  oh,  -  co^ 
=  Q(t)  . 


1 


(4.9) 


Let  x(t)  =  (Pi(t),  p2(t),  4>(t))  be  a  homoclinic  orbit  for  the  (p2,  p2,  (J))-dynamics  in  Pj  +  p2  +  p3 
=  M,  p3  =  J  where  M,  J  are  given  constants. 

As  in  Holmes  and  Marsden  [1983],  if  we  set 

y(t)  =  f  Q(s)  ds  +  y0  , 

Jo 
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we  have  only  to  prove  that  the  Melnikov  function 

M(v0)  =  J  {h0,  §-}  [x(t),  £  Q(s)  ds  +  Vo>  j)  dt  (4.10) 

has  simple  zeroes  to  get  "horseshoes",  where  {  ,  }  is  the  bracket  in  the  varibles  (jij,  1 x2,  (J>)  in 
(Hi  +  H2  +  li3  =  M>  M-3  =  J)- 

Note  that  £2  is  an  explicit  function  depending  on  |ii>  \i2,  and  so  will  not  be  a  constant  as 
time  changes.  For  reasons  that  will  become  clear,  we  will  consider,  instead  of  the  function  M(Vo), 
the  function 


i.e.. 


N(v0.  m>)  :=  M(v0)  *  yj-fe  ; 

'  L  { H°'  (x<0, 1 n(s)  ds + V0,  J) dt 


(4.11) 


Now, 


Hi'  * 


H, 


W3/I3 


(a(!2|ii  -  eabji2  cos  <J>)  cos(<J)  +  \\r)  +  b  cos  y(Ii|i2  -  mabe  cos  0)} 


From  now  on,  we  will  drop  on  Ij  and  I2,  remembering  Ij  ^  ea2  and  12  >  eb2. 
Let  us  first  consider  N(Vo« 0).  Then  £2(t)  =  -  (Ijp.2  -  eab  cos  <j>  |I])/A,  and 


¥-{ 


b  cos  v  +  a  cos(<j>  +  Vjt) 


I2Hi  -  eab(i2  cos  <j) ' 
IjP2  -  eabm  cos  (j> 


=  |  b  ^cos  J  £2(s)  ds  •  cos  Vo  ~  sin  J  £2(s)  ds  •  sin  Vo 


+  a 


I2P!  -  eab  JJ.J  cos  v) 
I2ji2  ~  Eabpj  cos  \j/ 


s(o  +  J[  fi(s) 

sin  ^<}H-  £  £2(s)  ds 


ds  cos  Vo 


sin  v 


(4.12) 


We  can  assume  without  loss  of  generality  that  Hi  and  [i2  are  even  functions  and  <j>  is  an  odd 
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function;  then  Q  is  an  even  function  and  so 


£ 


Q(s)  ds 


is  an  odd  function  of  t.  Thus 


( 


Hi' 

Ho'TT 


|  H0,  b  cos  J  Q(s) 


I2M-i  -  eab(i2  cos  <)> 
ds  +  a  t— - — — - r  cos 


Ijp.2  ~  £ab|i!  cos  <j> 


r  f‘ 

<j>  +  I  Q(s)  ds 
^  Jo  > 


{C1  -  eab|i2  cos  0 

H0,  b  sin  J  Q(s)  ds  +  a  7—: - ~rr . sin 


f  TT  I2H1  ~  eab}j.2  cos  <)> 

=  H0,a  — - — r- - — r  cos 


Ijli.2  -  eabjXj  cos  <j> 
I#!  “  £ab|X2  cos  <}> 


Ij^-2  -  eabiij  cos  <J> 

(J)  +  J  fi(s)  ds  j  j  cos  Vo 


£2(s)  ds 


[  H°’  a  I1|x2  -  eabji-!  cos  <]>  sin 


inW 


£2(s)  ds  f  sin  Vo 


By  symmetry,  the  first  term  of  (4.13)  will  vanish  after  integration.  Thus, 


N(\|/0,0)  =  -J  |h0. 


I2fij  -  eabji2  cos  (J) 
Ij^2  -  eabui  cos  ()>  SU1 


(♦+£ 


ft(s)  ds  f  dt  sin 


Thus,  what  we  have  to  do  is  prove  that 


£ka 


I2M-i  -  eab(i2  cos  $ 
Iijj-2  -  eabfij  cos  <)>  Sm 


4>+  f  fl(s) ds 

Jo 


dt  st  0 


Assume,  for  example,  that  Ij  =  I2  =  I  and  a  =  b  =  1.  Then  the  integrand  becomes 


tt  Ipi-e^coscJ)  .  (  ,  ^ 

Ho.  Tn.-c.^c^  sin  k  +  J0  n(s) ds 


Ip2  -  £|J-i  cos  (J) 

.  a  [  Im  -  £|I2  cos  (j)  . 
~  ^  9<}>  \  I|i2  -  ep.!  cos  $  sin 


♦+r 

,  Jo 


£2(s)  ds 


cos  Vo 
\  sin  Vo 


(4.13) 


(4.14) 


(4.15) 


and 
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Here 


9  el  sin  0  (H2  -  Hi )  . 

3$  {  >  "  .  . . ^  sm 


(IjJ-2  -  ejj-!  COS  <1>) 


-  ep2  cos  0 

+  - —  COS 

lM-2  -  EM-1  COS  0 


0  +  J  Q(s)  ds^ 


^0+J  Q(s)  dsj 


(4.16) 


(&-&)  >  =  “ 


(I2  -  e2  cos20)(^!  +  |i2) 


(I|i2  -  ejj-i  cos  0) 


sin  0  +  Q(s)  ds 


(4.17) 


Let  us  express  |il5  |I2  as  functions  of  0  for  one  of  the  homoclinic  orbits  which  is  in  the 
(Hj,  |j.2, 0)-dynamics.  It  is  given  as  the  intersection  of  (ij  +  (^  =  M  -  J  and 


„  1  2  T  2  „  .  v  1  (M-jr 

H°  =  2I (I^  +  “2e«*  0  II1II2)  =  4  . . j-"e . 


(4.18) 


(the  energy  of  the  homoclinic  orbit). 

Set  v  =  |x2  >  M-i>  Then  we  can  write 


H0 


4 


'  (M-J)2 
V I  +  e  cos  0 


+ 


I  -  e  cos  0 . 


(4.19) 


When 


H0  = 


1  f  (M-J)2 
4  l  i-e 


we  get 


v  =  ±  Vr  (M-  J)  where  T  =  —  _  (1  +  cos  0) 


1  +  f}  COS  0 


and  p 


Consider  the  homoclinic  orbit  obtained  from  v  =  V  T  (M  -  J).  Then 
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M-2  = 


i+Vr 


(M  -  J)  and  (ij  = 


\-<f 


(M-J)  . 


(4.21) 


Substitute  (4.21)  into  (4.16)  and  (4.17)  to  get 

9  _  _ 4  ft  sin  4>  Vr _ 

^  {(1  -  P  cos  <}>)  +  VF  (1  +  P  cos  <{>)}2 

(1  -  P  cos  0)  -  V  r  (1  +  P  cos  <{>) 
(1  -  P  cos  0)  +  V  T  (1  +  P  cos  <j>) 

4(1  -  p2  cos2(j))(M  -  J)"1 


in  +J  Q(s)  ds^j 

:°s  (♦+_[ 


Q(s)  ds  , 


f — _ {  t  }  _  _  _ 

v^2  ^1/  {(i  -  p  cos  <}>) + Vr  (l  +  p  cos  4>)}2 


sin 


(-1 


£2(s)  ds 


(4.22) 


(4.23) 


And 


a  1  tn  r  A 

t1!  = 


Thus  N(\}t0,  0)  =  Ax  +  A2  +  A3,  where 


.  -  r 


-  4P  sin 


((1  -  P  cos  <)>)  +  VF  (1  +  P  cos  <t>))2 


sinH 


£Xs)  ds 


dt  , 


=  _  J~  d-P cos  4>)  -  ^  r  a  +  p  cos  40  .  _  f  A ,  fl 


<j)  COS  I  4>  +  f  £Xs) 
(1  -  p  COS  0)  +  V  r  (1  +  p  COS  <())  V  Jo 


ds  dt  , 


(4.24a) 


(4.24b) 


w: 


1  -  P2  cos2<)> 


(d  -  p cos  <t>)  +  V r  (i  +  p  cos  <t>))2 


-j=  <j)  sin  ^<J)  4-  J*  £Xs)  ds  j  dt  .  (4.24c) 


Now,  note  that 


Q  =  -  —  (I|i2  -  e  cos  <J>  Hj)  = 


lM-2  ~  £  cos  4>  m 
I2  -  e2  cos2<J) 


1_  \h  -  P  cos  <b  lit 

*  1  -  p2  COS2<J) 


(4.25) 


and  that  all  the  integrand  decays  exponentially  as  t  — >  ±  «>.  Therefore  all  the  above  integrals  are 
analytic  functions  with  respect  to  l/l  and  P,  0  <  p  <  1,  0  <  1/1,  and  are  continuous  in  the  range 
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0  <  p  <  1,  0  <1/1.  We  will  consider  the  limiting  case  when  I  =  «».  Then  £2^0  and  so  (4.24) 
becomes 


.  -  /: 


-  4P  sin  V  T  <j> 


((1  -  P  cos  <j>)  +  VF  (1  +  P  cos  <t>))2 


sin<j>dt  , 


-[ 


(1  -  P  cos  <J>)  -  Vr  (1  +  p  cos  <j))  . 


s-L 


(1  -  P  cos  <}>)  +  Vr  (1  +  P  cos  <J>) 

1  —  p2  cos2<}> 

((1  —  P  cos  <J>)  +  Vr  (1  +  P  cos  <()))2  ^ 


<j>  cos  dt  , 


r  . 

=  <()  sin  0  dt 


Changing  variables,  these  become 


J-7C 


-  4p  sin20  Vr  d<p 


((1  -  P  cos  <|>)  +  V  r  (1  +  P  cos  <j>))2 
n  (1  -  p  cos  0)  -  Vr  (1  +  P  COS  (j>) 


a3 


_  r 

J-  Jt  (1  -  P  cos  <t>)  +  Vr  (1  +  P  COS  0) 

f  1  -  P2  cos2<t> 

J- Jt 


cos  <()  d(J)  , 


((1  -  P  cos  <t>)  +  Vr  (1  +  P  cos  (J)))2 


r 

sin  0  d()> 


(4.26a) 


(4.26) 


(4.26c) 


(4.27a) 


(4.27b) 


(4.27c) 


It  appears  to  be  difficult  to  directly  check  whether  the  sum  Aj  +  A2  +  A3  is  not  equal  to 
zero.  To  deal  with  this,  let  us  compare  the  order  of  Ax,  A2,  and  A3  as  V  P  — >  0.  The  order  of 
A!  is  0(p3/2).  For  A2  and  A3,  set  V  P  =  5  to  get 

(1  -  82  cos  <J>) - =  +  cos  <J»)(1  -  84  cos2<j>) 

Vl-52 

A2  =  -  I  - r .  . COS  <|)  dtp  .  (4.28) 

“  n  (1  -  52  cos  <J>)  +  —  ■  — . -  v  (1  +  cos  <t>)(l  -  84  cos2<J>) 

sn? 
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-L 


1  _  8  /  (1  + cos  $)(*  + 

*  l  -  82  V  1-82co 


O  COS  <{)) 
COS 


1  + 


V 


_8 _  /  (1  cos  4>)(1  + 

1-S2c< 


cos  <|>  d(J> 


8Z  cos  <1> 
cos  <J> 


The  integrand  is  equal  to 


COS<j) 


1  _  5  ^1  +  1  52  +  o(84)^(1  +  82  cos  <)>  +  0(84)Wl  +  cos<|> 

1  +  g  +  I  52  +  o(84))(l  +  82  cos  <j>  +  0(S4)W  1  +  cos  <j> 

=  (l  -  (8  +  j  83  (1  +  2  cos  <}>)  +  0(85))vi  +  cos<|>) 

x  { (l  -  8  (l  +  j  82(1  +  2  cos  <>)  +  0(84)^V  1  +  cos  0 J 

2 

+  -j  (1  +  cos  (())  82  ^1  +  ^  82(1  +  2  cos  <]>)  +  0(84)^ 

+  |  (1  +  cos  <(>)3/2  83  (l  +  ^  82  (1  +  2  cos  <j>)  +  0(S4))  |  cos  <|>  +  0(84) 
<}>  1 1  -  2V 1  +  cos  <|>  8  +  y  (1  +  cos  <]))82  j  +  0(83)  . 


=  cos 


Now  check  the  coefficients  of  8  and  82  in  (4.29);  the  coefficient  of  8  is  equal  to 


f71  r——~~  8 

2J  V 1  +  cos  <)>  cos  <J>  d<>  =  • 


For  A3, 


r  = 


v  7 

|3  -sin<j>(l-J3  cos  <J>)  -t-  2(3  sin  4>  cos  <}>  (1  +cos<)>) 
1  ~P  (1  +  Pcos<[>)2 


(4.29) 


(4.30) 


Thus  the  coefficient  of  8=V  p  in  the  expansion  of  A3  is  computed  as 
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-  sin2<J> 

V 1  +  cos  <f> 


Unfortunately,  the  first  order  term  in  8  in  (4.29)  vanishes,  so  we  must  check  the  coefficient 
of  S2.  From  (4.30),  the  coefficient  of  8s  in  A2  equals 


3  f* 

Y  J  cos  <{>  (1  +  cos  <J>)  d<{>  =  —  6tc  . 


Similarly,  the  coefficient  of  S2  in  the  expansion  of  A3  is  computed  to  be 


f*  2 

I  2  sin  <j)  d<}>  =  4k  . 
J-n 


Thus  the  coefficient  of  82  in  Aj  +  A2  +  A3  is  given  by  -  67t  +  47t  =  -  2te  *  0.  Therefore  the 
coefficient  of  52  in  the  expansion  of  N(y0>  0)  with  respect  to  P  when  I  — >  is  not  equal  to 
zero.  From  this,  we  can  conclude  that  the  Melnikov  integral  has  only  simple  zeros  for  generic 
parameter  values  if  the  distance  between  the  hinge  points  and  the  center  of  mass  of  the  third  body  is 
small.  We  summarize: 


Theorem  4.1  If  the  distance  between  the  center  of  mass  of  the  third  body  and  its  hinge  point  is 
sufficiently  small,  then  apart  from  isolated  values  of  the  system  parameters,  the  dynamics  of  the 
three  body  system  has  Poincari-Birkhoff-Smale  horseshoes,  so  is  non-integrable. 


§5  Discussion 


In  this  paper  we  have  developed  a  fairly  complete  picture  of  the  dynamics  of  the  planar  3 
body  stystem.  The  kind  of  analysis  we  have  presented  enables  one  to  study  equilibria  and  their 
stability,  bifurcations  of  equilibria,  periodic  and  chaotic  solutions.  Computer  graphics  of  the 
dynamics  illustrating  these  features  has  been  developed  by  Sreenath  [1987]. 

While  this  analysis  may  be  difficult  to  extend  to  a  complex  structure  of  n  bodies,  the  detailed 
understanding  of  the  dynamics  of  3  bodies  helps  us  to  understand  the  relation  between  chaos, 
coherence,  and  stability  in  more  complex  structures  and  in  the  continuum  limit  n  -»  In  fact,  it  is 
proved  in  a  preprint  of  Y.-  G.  Oh  that  the  straight  -  out  configuration  of  the  finite  coupled  rigid  n  - 
body  system  is  always  stable,  but  its  continuum  analogue  turns  out  to  not  be  formally  stable. 
Moreover,  we  also  have  a  good  understanding  of  the  special  structure  of  periodic  orbits  bifurcated 
from  the  equilibrium  for  the  symmetric  coupled  rigid  body  system. 

We  also  point  out  that  the  detailed  understanding  of  the  Hamiltonian  structure  via  symmetry 
reduction  should  assist  in  the  development  of  numerical  algorithms  and  the  control  theory  for  these 
systems. 

In  Grossman,  Krishnaprasad,  and  Marsden  [1987],  the  dynamics  of  coupled  three 
dimensional  rigid  bodies  is  studied.  The  analysis  there  indicates  that  there  may  be  a  symmetric 
Hamiltonian  Hopf  bifurcation  leading  to  interesting  periodic  and  chaotic  motions.  Again  one  can 
conjecture  the  possibility  of  interesting  three  dimensional  waves,  such  as  helical  waves,  being  built 
from  an  understanding  of  the  few  degrees  of  freedom  situation.  An  eventual  goal  is  to  link  this 
theory  up  with  the  infinite  dimensional  case  in  Krishnaprasad  and  Marsden  [1987],  Simo, 
Krishnaprasad,  and  Marsden  [1987]  and  Krishnaprasad,  Marsden,  Posburgh,  and  Simo  [1988]. 
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